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ABSTRACT 

Young stars are expected to gain most of their mass by accretion from a disk that forms around 
them as a result of angular momentum conservation in the collapsing protostellar cloud. Accretion 
initially proceeds at high rates of 10~ 6 — 10~ 5 M Q yr _1 resulting in strong irradiation of the stellar 
surface by the hot inner portion of the disk and leading to the suppression of the intrinsic stellar 
luminosity. Here we investigate how this luminosity suppression affects evolution of the protostellar 
properties. Using simple model based on the energy balance of accreting star we demonstrate that 
disk irradiation causes only a slight increase of the protostellar radius, at the level of several per cent. 
Such a weak effect is explained by a minor role played by the intrinsic stellar luminosity (at the time 
when it is significantly altered by irradiation) in the protostellar energy budget compared to the stellar 
deuterium burning luminosity and the inflow of the gravitational potential energy brought in by the 
freshly accreted material. Our results justify the neglect of irradiation effects in previous studies of 
the protostellar growth via disk accretion. Evolution of some other actively accreting objects such as 
young brown dwarfs and planets should also be only weakly sensitive to the effects of disk irradiation. 

Subject headings: stars: formation - accretion, accretion disks - planets and satellites: formation - 
solar system: formation 



1. INTRODUCTION. 

It is currently established that circumstellar disks are 
quite ubiquitous around young stellar objects of all 
masses (Muzerolle et al. 2003; Cesaroni et al. 2007). 
They represent an important ingredient of the star for- 
mation since initially protostars must be growing pre- 
dominantly by accretion through the disk: angular mo- 
mentum conservation forces the infalling protostellar 
cloud material to form a centrifugally supported disk 
which then accretes onto a star. It is quite likely that 
only a small fraction of the final stellar mass gets ac- 
quired by the direct infall onto the protostellar surface, so 
that almost all of the stellar mass gets processed through 
the disk. 

This picture of star formation implies very high initial 
mass accretion rates in the disk, at the level of 10~ 6 — 
10~ 5 M Q yr _1 as the Solar- type stars are thought to gain 
most of their mass during the first several 10 5 yr. The 
presence of such a high-M accretion flow just outside the 
protostar immediately raises an issue of its possible effect 
on the protostellar properties. 

There are several ways in which disk accretion affects 
the protostar. First, star gains mass from the disk which 
increases its binding energy and tends to make the star 
more compact. Second, accreting gas brings in some 
amount of thermal energy with it which contributes to 
the pressure support in the star. The exact amount of 
heat advected into the star with the accreted material 
is unknown but it seems likely that because of the disk 
geometry the accreted gas would have enough time to 
radiate away most of its thermal energy and would join 
the convective interior of the star with temperature much 
smaller than the stellar virial temperature. 

Third, intense energy dissipation taking place in the in- 
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nermost parts of the accretion disk leads to strong irradi- 
ation of the stellar surface by the disk (Frank & Shu 1986; 
Popham 1997). It has been recently realized (Rafikov 
2007) that irradiation can act to suppress the internal 
luminosity of the protostar similar to the suppression of 
the cooling of hot Jupiters by the radiation of their par- 
ent stars (Guillot et al. 1996; Burrows et al. 2000; Baraffe 
et al. 2003; Chabrier et al. 2004). Disk accretion at rates 
M ~ 10~ 6 — 10~ 5 M yr -1 can easily reduce internal 
stellar luminosity by a factor of several which may have 
important implications for the early stellar evolution. 

The goal of the present paper is to assess these impli- 
cations by following the evolution of accreting protostar 
properly taking into account effects of disk irradiation. 

2. METHOD OF CALCULATION. 

We consider a protostar of mass M and radius R grow- 
ing by accretion of gas at rate M from the circumstellar 
disk. In this work M is specified as an explicit function 
of time so that M(t) is also known. We assume the disk 
to extend all the way to the stellar surface as even a 
1 kG magnetic field (typical value measured in the ma- 
ture T Tauri systems, see Bouvier et al. 2007) would not 
be able to truncate the disk accreting at a high rate of 
M - 10~ 6 - 10~ 5 M yr -1 . At the stellar surface accret- 
ing gas passes through the boundary layer (Popham et 
al. 1993) in which its speed is reduced from the Keplerian 
velocity to the velocity of the stellar surface. 

To evaluate the effect of disk irradiation on the proto- 
stellar evolution we use an approach based on the energy 
conservation which was developed in Hartmann et al. 
(1997). In this approach the convective part of the star 
comprising most of its mass is assumed to behave like 
a polytrope with index n = 3/2, so that inside the star 
pressure P is related to the density p via P oc p 5 / 3 . This 
approximation works very well in highly ionized, dense, 
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and fully convective interiors of young stars. The total 
energy of such a star (a sum of its thermal and gravita- 
tional energies 2 ) is E to t — — (3/7)GM 2 /i? (Kippenhahn 
& Weigert 1994). Evolution of protostellar properties 
- luminosity L, radius R - as a function of time [or, 
equivalently, stellar mass M(t)] is then governed by the 
following equation: 
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The l.h.s. of this equation represents the change in the 
total stellar energy, the first and second terms in the 
r.h.s. are the gravitational potential energy and the ther- 
mal energy brought in with the accreted material, while 
Ld(M, R) is a deuterium luminosity of a protostar. Stel- 
lar luminosity L is the luminosity carried towards the 
photosphere by the convective motions in the stellar in- 
terior. It is different from the integrated emissivity of 
the stellar surface since the star also intercepts and rera- 
diates a fraction of energy released in the accretion disk. 
Rate at which thermal energy gets accreted by the star 

is 
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where 7 is the ratio of specific heats of accreted gas 
(which can be different from 7 = 5/3 characteristic for 
the stellar interior), k B is the Boltzmann constant, \x and 
T are the mean molecular weight and the temperature 
of the accreted gas. Dimensionless parameter a can be 
written as 
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where T V i r is the stellar virial temperature. Gas accret- 
ing from the disk experiences strong dissipation in the 
boundary layer near the stellar surface. In this layer 
gas temperature can become an appreciable fraction of 
T v i r . However, the cooling time in the boundary layer 
and the outermost layers of the star is very short so that 
the accreted gas cools efficiently and should ultimately 
join stellar interior with temperature T which is much 
lower than T V i r [unless M is extremely high, in excess 
of 1CT 4 M Q yr" 1 , see Popham (1997)]. Thus, under the 
conditions considered in this work one expect a <C 1 
but the actual value of this parameter is rather poorly 
constrained (it depends on the gas thermodynamics, ra- 
diative transport, and the viscosity prescription in the 
boundary layer which are poorly constrained). As here 
we are primarily interested in the effects of disk irradia- 
tion we set a = for simplicity 3 . In this case equation 
(JTJ can be rewritten as 
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This is an evolution equation for R and can be easily 
integrated numerically once the dependencies of Ld and 
L on stellar parameters are known. 

2 In this work we neglect stellar rotation. 

3 Prialnik & Livio (1985) and Hartmann et al. (1997) have previ- 
ously investigated the effect of the variation of a on the protostellar 
evolution. 



For Ld we adopt the expression obtained in Stahler 
(1998) by integrating the rate of energy release due to D 
burning within the n = 3/2 polytropc: 
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where Ld,o = 1-92 x 10 17 L Q and fn is the fractional 
D abundance relative to the initial D number abundance 
[D/H] taken to be 2 x 10 -5 . Parameter fu is not constant 
in time - it evolves since D burns in the stellar interior 
while the new D is being brought in with the accreting 
material (with the initial abundance [D/H]). As a result, 
one finds (Stahler 1988; Hartmann et al. 1997) 
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where (3d — 9.2 x 10 13 ergs g 1 is the energy released by 
D fusion per gram of stellar material (assuming [D/H] = 
2 x 1CT 5 ). 

The most important aspect of this work which distin- 
guishes it from Hartmann et al. (1997) is the calculation 
of L. In the approximation adopted by Hartmann et al. 
(1997) L is a function of R and M only. In our case sit- 
uation is different: irradiation of the stellar surface gives 
rise to an outer convectively stable layer below the stel- 
lar photosphere (Rafikov 2007), similar to the radiative 
layer that form in the atmospheres of the close-in giant 
planets irradiated by their parent stars (Guillot et al. 
1996; Burrows et al. 2000). This external radiative zone 
suppresses the local radiative flux coming from stellar in- 
terior and this changes the integrated stellar luminosity, 
which becomes a function of irradiation intensity. As a 
result, in irradiated case L depends not only on R and 
M but also on M (which determines the strength of the 
irradiation flux). 

Rafikov (2007) has demonstrated that for a given opac- 
ity behavior at the stellar surface (parametrized in his 
case to be a power-law function of gas pressure P and 
temperature T, n oc P a T^) the degree of luminosity sup- 
pression depends only on the so-called irradiation param- 
eter 
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which is (up to a constant factor) the ratio of the accre- 
tion luminosity of the disk GMM /R to the luminosity 
Lq that a star would have had in the absence of irra- 
diation. Suppression factor x(A) -4 1 as A <C 1 while 
x(A) < 1 when A > 1. In a simple case considered 
by Rafikov (2007) the dependence of \ on the opacity 
behavior comes only through the parameter 4 



i = P+{l + a)/V ad , 
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where S/ a d is the adiabatic temperature gradient near the 
stellar surface. 

In the strongly irradiated case the temperature of the 
stellar surface varies as a function of latitude: equatorial 
belt is strongly heated by the hot inner parts of the disk 
while the polar regions of the star are virtually unaffected 

4 For convection to set in at some depth below the outer radiative 
zone a condition £ > 4 must be fulfilled, see Rafikov (2007) for 
details. 
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by irradiation and preserve their temperature at the level 
of To = (Lo/4vri? 2 cr) 1 / 4 . In this situation one may won- 
der whether it is reasonable to assume a fixed opacity 
law (as was done in Rafikov 2007) for the calculation of 
luminosity suppression given that the behavior of n can 
be different between the polar and the equatorial regions 
of the star. However, is was shown in Rafikov (2007) 
that for the opacity behavior typical for stellar photo- 
spheres in the temperature interval from ~ 2.5 x 10 3 K 
to 10 K the cooling of irradiated stars occurs mainly 
through their polar caps 5 (even though the opacity scal- 
ing with temperature and pressure changes quite drasti- 
cally within this temperature interval at around 5000 K). 
As a result, no matter how hot the equatorial parts of the 
star become and how complicated the opacity behavior is 
in this portion of the stellar surface, the integrated stel- 
lar luminosity L does not depend on these details very 
strongly but is rather determined by the properties of the 
polar regions of the star: the size of the cool polar caps 
in which the photospheric temperature is preserved at 
the level of To, and the opacity behavior in the adjacent 
parts of the stellar surface. This provides motivation for 
using the stellar luminosity prescription represented by 
equation |7|. In this paper we adopt k characteristic for 
the temperature interval 2.5 x 10 3 K < T < 5 x 10 3 K 
(Bell & Lin 1994): 

k « 6 x io-i 4 p 2 / 3 T 7 / 3 . (9) 

This scaling should be reasonable in the polar regions 
of irradiated stars where the photospheric temperature 
To is not strongly affected by irradiation and is close 
to the photospheric temperature that an isolated non- 
accreting star would have possessed in the Hayashi phase. 
The opacity law (J9j> corresponds to £ = 6.5, assuming 
Vad = 2/5 as appropriate for the fully ionized stellar 
interior (7 = 5/3). Having specified the opacity law we 
have thus fully determined the behavior of x(A) (which 
we take from Rafikov [2007] , see the curve corresponding 
to £ = 6.5 in Fig. 4 of that paper) necessary to compute 
L. 

One remaining ingredient of the calculation is the 
choice of Lo(M, R) - the luminosity of a non-irradiated 
star. Hartmann et al. (1997) have adopted the following 
fit to the stellar evolution tracks of D'Antona & Mazz- 
itelli (1994): 

'W^- 1 Mans) Km) ' (10) 

Evolution tracks in DAntona & Mazzitelli (1994) have 
been calculated using the equation of state from Magni 
& Mazzitelli (1979) that has been superceded by the 
more refined treatments (Saumon et al. 1995). Also, 
DAntona & Mazzitelli's treatment of convection is based 
on Canuto & Mazzitelli (1991) which has previously 
raised some concerns (Demarque et al. 1999; Nordlund & 
Stein 1999). Despite these deficiencies, we have chosen 
to adopt the prescription (fl"0|) in our work because of its 
simplicity 6 and also to allow direct comparison with the 
results of Hartmann et al. (1997). 

5 Situation is different in the case of giant planets irradiated by 
the circumplanetary disks, see Rafikov (2007) for details. 

6 We do not require perfect knowledge of Lq since our primary 
goal is to evaluate the importance of disk irradiation. 




M [M s ] 

Fig. 1. — Evolution of the protostellar properties for M = 
4 X 10~ 6 M Q yr" 1 . Initially M = 0.1 M and R = 1.5 R . 
We display the runs of irradiation parameter A (a), suppression 
factor x(A) (b), ratio Lj /Lq of the D burning luminosity to the 
luminosity given by Eq. q 10 ^ (c), stellar radius R (d), and the D 
abundance / rj relative to its initial value ( e) , as functions of stellar 
mass M. In the last three panels we display the corresponding 
quantities both in irradiated (solid) and non-irradiated (dashed) 
cases. 



Equations (0, (©, ©, ©, & |10|) supplemented with 
dM/dt — M(t) and the dependence x(^) from Rafikov 
(2007) fully determine the evolution of an accreting pro- 
tostar irradiated by its own disk. This system of equa- 
tions is then evolved numerically assuming that the pre- 
scription for M{t) is given. 

3. RESULTS. 

Here we present the results of our calculations. As 
initial conditions we choose M = 0.1 M© and fn = 1. 
We vary R and the prescription for M(t) to see their 
effect on the evolution of stellar properties. 

In Figure [1] we display protostellar evolution for ini- 
tial R = 1.5 R Q and a uniform M = 4x 10~ 6 M Q yr -1 
with and without the effects of disk irradiation included. 
One can see that in both irradiated and non-irradiated 
cases evolution is pretty much the same: star initially 
contracts until its central density and temperature be- 
come high enough for the deuterium to ignite. Right 
after that D burning strongly dominates over the stellar 
luminosity and D abundance /d starts going down. Re- 
sulting energy release in the stellar interior causes star 
to expand out to 2.5 R© and D luminosity decreases 
appreciably (but still exceeds Lq by a factor of several). 
These results are in full agreement with the calculations 
of Hartmann et al. (1997). 

In the top two panels of Figure [1] we show quanti- 
ties unique for the irradiated case: run of the irradia- 
tion parameter A and the suppression factor x(^). As 
R initially decreases A goes up to « 1.5 x 10 3 since 
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Fig. 2. — Relative difference between the stellar radii in the 
irradiated R and non-irradiated Rq cases (panel c) for a protostar 
accreting at a uniform rate M = 2 X 10 -6 Mq yr — x . Different 
curves correspond to different initial radii: 1.0 R.0 (solid), 1.5 Kq 
(dashed), and 2.0 R© (dotted). Panels a and b display runs of 
irradiation parameter A and suppression factor \ as a function of 
stellar mass. 



A oc MM 01 R~ 3 - 34 for L given by equation (JTOj). As a 
result, the internal stellar luminosity is appreciably sup- 
pressed and x reaches ks 0.55 demonstrating the impor- 
tance of disk irradiation in regulating L. 

At the same time, although the effect of irradiation on 
L is of order unity, Figure Q] clearly demonstrates that 
irradiation affects other stellar properties such as R and 
fo only very weakly. As expected, R in irradiated case 
is larger than in the non-irradiated case (because L sup- 
pression allows more heat to be retained inside the star) 
but only by several per cent. 

To check that this result is not an artefact of our initial 
conditions and assumed M we have additionally calcu- 
lated protostellar evolution for uniform M = 2 x 10~ 6 
M Q yr^ 1 and M = 10~ 5 M Q yr -1 and different initial 
R. Results are presented in Figures [2] and [3] in which we 
display the evolution of AR/R — (R — R n )/ R , where R 
and i? are the values of the protostellar radius with and 
without disk irradiation taken into account. It is quite 
clear from these plots that despite the rather severe lu- 
minosity suppression (x reaching 0.3 in the high M case 
for the initial R = 1 R©) the relative stellar radius in- 
crease due to irradiation is always rather small. Note 
that at M ~ M Q we find AR/R « 4% when M = 10~ 5 
M© yr -1 (Figure [2]) which is larger than in the higher M 
case shown in Figure [3] (when AR/R ks 3%). 

Very similar outcome has been found in the case of 
non- uniform M: in Figure |4] we plot AR/R for different 
initial R and M — aM with a chosen in such a way 
that a protostar grows to 1 M© in I0 5 yr. This growth 
time is close to the time needed to reach 1 M© in the 
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Fig. 3.— Same as Fig|2]but for a uniform M = 10~ 5 M Q yr -1 . 



case of constant M = 10 5 M yr 1 , see Figure [3J For 
non-uniform M one again finds that irradiated protostar 
differs in radius from the non-irradiated star by only a 
few per cent (AR/R w 3% at M ~ M Q ). It is also 
obvious from the comparison of Figures [3fc and 0J; that 
roughly the same growth time translates into very similar 
behavior of AR/R in the two cases independent of how 
M(t) evolves. We thus conclude that irradiation by the 
circumstellar accretion disk has rather small effect on 
the protostellar evolution and that this outcome is rather 
generic. 

4. DISCUSSION. 

Weak sensitivity of R to disk irradiation given the 
strong effect that irradiation has on the internal lumi- 
nosity L of a protostar may seem surprising. However, 
one should bear in mind that besides L there are other 
contributors to the stellar energy budget, namely the D 
burning luminosity Lrj and the gravitational potential 
energy gained with the accreted material. It turns out 
that these contributions dominate the energy budget over 
L. 

To see this we rewrite equation Q using definition J7]) 
in the following form: 
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In the right-hand side of this equation the first term in 
brackets describes the relative role of D burning in the 
total energy budget, second term represents intrinsic stel- 
lar luminosity, while the third term is due to the inflow 
of the gravitational potential energy GMM/R - the ra- 
tio of this energy inflow to Lq differs from A only by a 
constant factor. 

Suppression of L is largest when A is highest which is 
easy to see by inspecting Figures [T]|4] But this automat- 
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Fig. 4. — Same as Fig [2] but for M oc M and the growth time 
10 5 yr to reach M = 1 Mq. 



ically means that the maximum deviation of the sup- 
pression factor x from unity occurs precisely when the 
inflow of the gravitational potential energy far exceeds 
the stellar luminosity. Apparently, under these circum- 
stances L is a subdominant contribution to the stellar 
energy budget and thus even a significant reduction of L 
compared to Lq is going to be negligible compared to the 
gravitational energy influx. 

Moreover, Figure Q] demonstrates that A reaches its 
maximum when R is at its minimum while fr> is still very 
close to unity. At this point vigorous D burning com- 
mences inside the star giving rise to very high Ljj/Lq. 
As a result, at the evolutionary stage when x is minimal 
L is subdominant in comparison to not only GMM/R 
but also L d . This additionally downplays the role of the 
luminosity suppression by irradiation in the early proto- 
stellar evolution. 

This line of reasoning also explains why at M ~ 1 M© 
we have found AR/R to be larger for lower M (see JJ]). 
First, smaller M means lower A so that the ratio of L to 
the gravitational energy inflow rate GMM/R in the low 
M case is larger than in the high M case. Also, at M ~ 
M Q one generally finds /b < 1 (see Figure QJ) so that 
Ld is mainly due to the burning of the freshly accreted D 
(rather than the D that remained in the protostar from 

previous accretion). Since in the low M case less fresh 
D is supplied to the protostar Ld must also be lower 
than in the high M case. As a result, in the lower M 
case L plays a more significant role compared to Lq (in 
which case AR/R should be more sensitive to changes in 
L caused by irradiation) than in the high M case. 

This conclusion immediately raises the following ques- 
tion: since AR/R increases as M decreases would one 
find AR/R ~ 1 at low enough Ml The answer is no, 
and it has to do with the fact that x appreciably differs 



Fig. 5. — Plot of A50 (the value of A at which the suppression 
factor x = 0-5) as a function of parameter £ characterizing the 
opacity law in the outer layers of irradiated protostar. Note that 
A50 never decreases below several hundred. 



from unity (obviously, a necessary condition for getting 
AR/R ~ f) only at rather high A. This is a generic 
feature of disk irradiation which is illustrated in Fig- 
ure [5] where we display A50 - the value of A at which 
x(A5o) = 0.5 - as a function of assumed opacity law rep- 
resented by the parameter £, see equation (jSJ). One can 
see that A50 > 10 2 for all £ > 4, meaning that significant 
luminosity suppression requires rather high M. This in- 
efficiency of irradiation in suppressing L is caused by the 
specific geometry of disk irradiation in which the irradia- 
tion flux is a very sensitive function (oc 9 5 ) of the latitude 
at the stellar surface 6, see Rafikov (2007). Because of 
that stellar polar caps can stay cool even at rather high 
M allowing unsuppressed flux to be emitted over a sig- 
nificant portion of the stellar surface. 

In the case of £ = 6.5 as appropriate for cool, low-mass 
protostars one finds that A50 = 2.2 x 10 3 which accord- 
ing to equation j7]) immediately implies that GMM/R w 
175-Lo when x = 0.5. Clearly, in this case stellar luminos- 
ity should have small effect on the protostellar evolution. 
If M is reduced 7 so that GMM/R ~ L Q (and A - I) 
stellar luminosity would be playing a significant role in 
the stellar energy budget, however x would be very close 
to unity (see Rafikov 2007) and the L suppression by 
irradiation would be negligible. Thus, under no circum- 
stances should one expect AR/R larger than several per 
cent, meaning that quite generally the irradiation by ac- 
cretion disk is unlikely to play a significant role in the 
evolution of the protostellar properties. 

Looking at this conclusion under slightly different an- 

7 At small enough M (M < 10 — 6 Mq yr —1 ) protostellar mag- 
netic field is likely to disrupt accretion flow (Konigl 1991) thus 
invalidating our assumption of the direct disk accretion onto the 
star. 
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gle, our results also imply that when considering evo- 
lution of the protostars accreting at M ~ 10~ 6 — 10~ 5 
Mq yr _1 one may completely neglect L in the stellar en- 
ergy budget and still get rather decent description of the 
protostellar evolution. 

Our results are obtained assuming a specific value of 
£ typical for the low-mass stars. However, one expects 
our major conclusions to remain valid also in the case 
of other accreting objects such as young brown dwarfs 
and giant planets. Although these objects are likely to 
be characterized by values of £ different from 6.5, our 
Figure [5] clearly demonstrates that the efficient suppres- 
sion of stellar luminosity requires very large A for any 
value of £. As a result, we generically expect L to be a 
subdominant energy contributor for any actively accret- 
ing object meaning that its radius (and other properties) 
should be different from the radius of its non-irradiated 
counterpart by at most a few per cent. 

Previous evolutionary studies of the protostellar as- 
sembly by accretion (Mercer-Smith et al. 1984; Palla & 
Stahler 1992; Siess & Forestini 1996; Hartmann et al. 
1997; Siess et al. 1997, 1999) have always neglected the 
effect of disk irradiation on stellar properties. Results 
of this work demonstrate that this simplification should 
not have led to major deviations from reality since disk 
irradiation affects stellar properties only at the few per 
cent level. Protostellar evolution may more sensitively 
depend 8 (Prialnik & Livio 1985; Hartmann et al. 1997) 
on the exact amount of thermal energy that gets accreted 
with the disk material which is parametrized by the pa- 
rameter a, see equation ([3]). Disk irradiation is likely 
to play a very important role is setting the value of a 

8 As long as a is not very small. Whether this is true is debatable 
(Popham 1997), and we assume a = in this work. 



since the conditions in the external radiative zone near 
the stellar surface regulate the thermodynamic proper- 
ties of the gas at the convective-radiative boundary of the 
protostar, which is likely to be crucial in determining a. 

5. SUMMARY. 

We considered evolution of the properties of a growing 
protostar accreting gas from the circumstellar disk with 
the goal of assessing the impact of irradiation by the in- 
ner regions of the disk. We find that disk irradiation 
plays minor role in the radius evolution of a protostar so 
that the radius of an irradiated star is larger by only a 
few per cent compared to the non-irradiated case. This 
result is largely independent of the specific behavior of 
the protostellar mass accretion rate as long as M is high 
enough, at the level of 1CT 6 - 10~ 5 M Q yr _1 . The weak 
sensitivity of the stellar properties to the disk irradiation 
is explained by the minor role played by the stellar lu- 
minosity L in the energy budget of the star at the time 
when L is significantly altered by irradiation. At these 
stages of protostellar evolution the inflow of the negative 
gravitational potential energy brought in by the accret- 
ing material and the deuterium burning luminosity are 
the much more significant contributors to the stellar en- 
ergy budget. Our results imply that the previous studies 
of the protostellar accretion which neglected disk irradi- 
ation may have underpredicted stellar radius only at the 
level of several per cent. Our major conclusions should 
be directly applicable to the evolution of other actively 
accreting objects such young brown dwarfs and planets. 
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